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ABSTRACT 

Requirements for magnetic coupling and accretion in the active layer of a pro- 
tostellar disk are re-examined, and some implications for thermal emission from 
the layer are discussed. The ionization and electrical conductivity are calculated 
following the general scheme of Ilgner and Nelson but with an updated UMIST 
database of chemical reactions and some improvements in the grain physics, and 
for the minimum-mass solar nebula rather than an alpha disk. The new lim- 
its on grain abundance are slightly more severe than theirs. Even for optimally 
sized grains, the layer should be at least marginally optically thin to its own 
thermal radiation, so that narrow, highly saturated emission lines of water and 
other molecular species would be expected if accretion is driven by turbulence 
and standard rates of ionization prevail. If the grain size distribution extends 
broadly from well below a micron to a millimeter or more, as suggested by ob- 
servations, then the layer may be so optically thin that its cooling is dominated 
by molecular emission. Even under such conditions, it is difficult to have ac- 
tive layers of more than lOgcm"^ near 1 AU unless dust is entirely eliminated or 
greatly enhanced ionization rates are assumed. Equipartition-strength magnetic 
fields are then required in these regions of the disk if observed accretion rates 
are driven by magnetorotational turbulence. Wind-driven accretion might allow 
weaker fields and less massive active layers but would not heat the layer as much 
as turbulence and therefore might not produce emission lines. 

Subject headings: accretion disks — magnetic fields — molecular processes — stars:pre- 
main sequence 



Introduction 



Obs ervations of classical T Tauri stars indicate accretion rates M ~ 10 ^^^MqJt ^ for 
lO^^Vr. fjHartmann et al.lll998l ). The torques driving the accr etion are probably magne tic, 



whether exerted on the face of the disk by a magnetized wind (IBlandford fc Paynelll982l ). or 
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across the thickness of the disk by magnetorotational (MRI) turbulence ( iBalbus fc Hawley 



199ll . Il998l ). At the low temperatures characteristic of most of the mass of a protostellar 



disk, the ionization fraction would be too low for good coupling to a magnetic field were it 
not for nonthermal ionization processes such as cosmic rays, X-rays, and radioactivity. Since 
the sources of the first two lie outside the disk, the surface layers of the disk are expected to 
be more strongly ionized than the midplane. It is possible that accretion occurs mainly in the 
surface layers, while the rest of t he disk column forms a magnetically decoupled "dead zone" 
in which mass may accumulate (jGammielll996l ). Dead zones may be favorabl e environments 



for the formatio n of planets, or for the survival of planets once formed (e.g. 
Ida fc LinI [20081 ). 
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It has been recognized at least since Gammie's work that small dust grains, through 
their influence on the ionization balance, are crucial for MRI turbulence and for magnetic 
coupling more generally. A number of authors have investigated the thickness of the active 
layer based on different disk rnodels, ion ization sources, ch e mical reaction networks, and grain 
properties. iGlassgold et al.l (119971 ) a ndllgea fc Glassgoldl (Il999l ) pointed out the importance 
of stellar X-rays to the active l ayer. ISano et al.l (l2000l ) considered the minimum mass solar 
nebula (MMSN, iHayashil Il98ll ) wit h a chemical n e twork including dust grains and found 
dead zones extending to ~ 15 AU. iFromang et al.l (120021 ) showed that the conductivity is 
sensitive to the abundance of "metals" (Mg, Fe, etc.) in the gas phase, eliminating the dead 
zone entirely from some o f their grain- fre e mod els, which were a disks with lower column 
density than th e MM SN. ISemenov et al.l (12004 ) adopted the "fiducial" a disk model by 
D'Alessio et al.l (|l999l ). which also has lower surface d ensity at 1 AU than the MMSN, and 
a chemical network based on th e UMIST 95 dat a base (IMillar et al.lll997l ) supplemented by 
reactions on grain surfaces fro mlHasegawa et al.l (119921 ): they found a marginally dead zone 
somewhat smaller than that of ISano et al.l (120001 ). More recently, in a similar study based on 
an a-disk model and a yet more extensive network, lllgner fc NelsonI ( l2006al ) concluded that 



adding sub-micron sized grains with concentration of = 10 per H2 molecule efficiently 
depletes metal atoms and dramatically reduces the extent of the active layer. 

Dust grains in protostellar disks have observable signatures in the spectral energy dis- 
tribution (SED) from infrared to millimeter wavebands. Small grains in particular near the 
disk surface are superheated by absorbtion of visible light from the star; this promotes flar- 
ing (concavity) of the disk surface and may contribute to th e relat ively slow decline of the 
SED with increasing wavelength (IChiang fc GoldreichI Il997l . Il999l ). More refined models, 
however, show that th e dust size distr i bution may not be well constrained solely from SEDs 
( IChiang et al.ll200ll ). iD'Alessio et al.l (120061 ) further studied comprehensively the effect of 
grain growth and settling on the SED of the protoplanet ary disks; comparing th eir theoreti- 



cal SED with recent mid- and near-IR observations (e.g., Hartmann et ahlbood ). they found 
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that grains in the disk upper layer should be depleted by at least a factor of 10 relative to 
interstellar dust-to-gas ratios, and perhaps by 10^ — 10'^. 

Little attention has been given to the implications of the conductivity constraints for the 
emissions by the dust and gas. This is the focus of the present paper. We recalculate, in yet 
greater detail, the influence of grains on the ionization balance. In view of the ambiguities in 
measuring small grains from the SEDs and the difficulty of improving on previous work on 
that problem, however, we chose to concentrate on an indirect signature of dust depletion: 
enhanced molecular emission from the active layer. The fundamental CO emission line has 
been frequen tly observed in CT TSs, with the likely origin of the protostellar disks within 



1-2 AU (e.g., iNajita et al.ll2003l ). Recently, discoveries of organic and water molecules have 



been reported, and these molecules are likely to originate froni the inner region (< 5 Au) of 
young protostellar disks jCarr fc Najitalbnosl : ISalvk et all 120081 : iMandell et allbnosl ). These 
molecular emissions are consistent with hot, optically thin gas. Based on the derived tem- 
perature and column density, the emission lines may originate from a UV-heated layer above 
the disk surface. However, we argue from the conductivity constraints that dust must be 
sufficiently depleted from the active layer so that it is at least marginally optically thin, 
allowing molecular emission lines to stand out. If the dust size distribution extends to grains 
as small as 10~^/im, then the heat of accretion may be expressed primarily in these lines 
rather than in the dust continuum. 

The organization of the paper is as follows. In §2], we review constraints on the minimum 
ionization fraction and magnetic field strength required to drive accretion at observed rates 
by winds or by MRI turbulence. In §§3-4, we calculate the thickness of the active layer. This 
is an extension of the work by IN06a, with the following differences and additions. Firstly, 
we use the MMSN for the density and temperature profile of the disk, rather than an a 
model. The MMSN is simply defined and has been used extensively as a standard model 
for discussion of protostellar disks and planet formation. It has an empirical basis in the 
observed properties of our own solar system and is not inconsistent with inferences of T Tauri 
disks by, e.g., submillimeter emission. Steady a disks are have a weaker observational basis. 
They come in many forms, depending on assumptions about accretion and cooling, and 
they are often not self-consistent where active layers exist. The particular a model used by 
IN06a is convex rather t han fiared like the M MSN. Secondly, we adopt the latest version 
of the UMIST database (jWoodall et al.l 120071 1 for the chemical reactions. Changes in the 
reaction rates have perceptible consequences for the ionization fraction. Thirdly, we provide 
more detailed analysis of the dependence of conductivity on the grain properties, especially 
their size distribution. In §5, we discuss the conditions under which molecular lines in the 
active layer of protostellar disks may stand out above the dust continuum. For illustrative 
purposes, we study the water molecule only; it is abundant, and its rich spectrum makes it 
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a good coolant. Using the emission line and energy level list given by iBarber et al.l (120061 ) 
(hereafter BT), we calculate the molecular line emissions from an isothermal slab in local 
thermodynamic equilibrium (LTE) to mimic the disk active layer. Various line broadening 
effects are considered in our calculation. We derive the conditions on the dust abundance 
and size distribution in which water lines dominate the cooling of the layer. Limitations of 
our study and some open questions are discussed in §6, and a summary is given in §7. 



2. Requirements for Magnetically-Driven Accretion 

In this section, we discuss theoretical constraints on the ionization fraction (xg), mag- 
netic field strength {Bz or -Brms), and surface mass density (Sa) of the active layer re- 
quired to explain magnetically driven accretion rates of order 10~^ M_7 yr~^, where 
M_7 = M/(lO-^M0yr-^). First we consider a magnetized wind, and then MRI turbu- 
lence. In the latter case, we emphasize the conditions needed to sustain accretion at the 
above rate, which are more stringent than those required for linear instability. The mini- 
mal field strength {B) is about one order of magnitude larger for MRI turbulence than for 
wind-driven accretion at the same M. But the degree of ionization is similar. 



2.1. Disk Model 



We adopt the minimum-mass s olar nebular (MMSN) model with stellar mass M* = IMq 
(jHayashilll981uHayashi et al.lll985l ). The disk is vertically isothermal with radial temperature 
profile 



surface mass density 
sound speed 

c 

volume density 



T = TorAu 
/ kT \ 



1/2 



To = 280K, 
So = 1700g cm-2 

'1/4. -1 



O.ggr'u' kms-S ^ 2.34, 



p{r, z) = po{r) exp{-z /2h ) , po{r) 



and scale height 



/iV27r 

h = Cs/^K = 0.03ri{f AU 



1.4 X 10 ^r^u^'^^gcm 



(2) 
(3) 

(4) 
(5) 
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where is the disk radius measured in astronomical units and Qk 
that the disk flares, that is, h/r increases with r. 



(G'M,/r3)V2. Note 



2.2. Magnetized wind 

When accretion is wind-driven, angular momentum is extracted via the mean magnetic 
torque per unit area —rBzB^/An exerted on the surface the disk. This must balance the 
rate of loss of angular momentum per unit area, which is —MQ/8tt when one considers that 
the accretion is divided between the two faces of the disk so that 



i\BMB^\)>Mn/2r. 



(6) 



We consider M > for inflow, and hereafter we omit the overbars and absolute value signs 
where this will not cause confusi on. It is also required tha t Br > B^/ in order that the 
wind be centrifugally propelled (IBlandford fc Payndll982l ). Minimizing the total magnetic 
ener gy subject to these constraints leads to a minimum total field at the surface of the disk 



[e.g. 



Wardlell2nn7h 



B > 



2 MQ 
V3 r 



1/2 



0.31 rlu^^ G 



(7) 



In th is configuration, Br = —B^/2; some wind models require \Br 
19931 ) and therefore a stronger total field for the same M. 



> \B^\ (e.g.. lWardle &: Koenigl 



The active layer has thickness ha ~ ^c^/fi, where the dimensionless factor ^ is such 
that the volume density at the base of the active layer is pa = I^a/ha- Since the vertical 
density profile under isothermal conditions is gaussian with scale height Cg/Q, C, depends 
logarithmically on E^/E; for example, if E^ = lOgcm"^, then the appropriate value of 
^ ~ 0.39 at 1 AU. The roughness of our estimates hardly justifies keeping track of this 
logarithmic dependence, so we simply take ha = O.Scs/O at all radii hereafter. 



In the configuration that achieves the minimum ([7]), the total horizontal component 
Bh = {B^ + B'^y^'^ is approximately equal to \Bz\ at the disk surface. But in the dead 
zone, Bh <^ Bz because the field there is poorly coupled to the gas and straightens out 
under magnetic tension. Thus in the active layer, \dBh/dz\ > \Bz\/ha 3> \dBz/dr\, so that 
there will be a horizontal elect ric curre n t in t he active layer and an associated electric field. 
As described in the review by IWardld (120071 ). the current and the electric field are not in 
the same direction under the likely conditions prevailing in active layers; the conductivity 
is nontrivially tensorial because of the field itself, and because some charged species are 
more firmly attached to the field than others. This is quantified by the Hall parameters 
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Pj = ojjtj for each charged species j, where Uj = \qjB\/mjC is the cyclotron frequency, and 
tj = TTijirnn/ Pjp{(yv) jn is the timescale on which particles of species j lose their momentum 
in collisions with the neutrals. The reduced mass pj = mjirtn/ {rrij + mn) is essentially me for 
free electrons, and pj ~ m„ 2mp for metallic or molecular ions as well as for charged grains, 
since all of these are heavier than molecular hydrogen. Wardle gives numerical scalings for 
free electrons and ions, which we rewrite in terms of Si = Ea/(10gcm~^), Bq = 5/(1 G), 
and the standard properties of the MMSN as 



/3, ^ 87 ^^'Bcr At , A ^ 0-19 ^^'Bgt 



AU 



Charged grains are unmagnetized, jSg <^ 1. For magnetic fields comparable to the value 
([7]), the ions follow the neutrals (because Pi <^ 1) but the electrons follow the field (because 
/3e ^ !)• In the rest frame of the neutrals, the current density J is therefore borne almost 
entirely by the electrons, as is the Lorentz force J X B/c. So there must be an electric field 
E ^ — J X B/eUeC in this frame to prevent the electrons from accelerating, since collisions 
are ineffective. Thus "Ohm's Law" in the frame of the neutrals becomes anE ^ J x 
where B = B/B is a unit vector and an ~ —eUec/B is the Hall conductivity. The electric 
field component parallel to B is much smaller than the perpendicular component because 
the ordinary Ohmic conductivity ~ e^riete/rrie = —Pe<^H is much larger than an- 

The electric field in the inertial reference frame of the star, where the neutrals have 
velocity v, is E' = E — v x B /c. The azimuthal component E'^ must vanish, at least on 
average, provided that there is no secular increase in the magnetic flux threading the disk. 
Substituting E ^ a^^^J X B ^ {B x V X B)/ATiene into E'^ = + {yrB^ - v^B^)/c, and 
presuming that f ^ ^ t",. <^ within the active layer (though not within the wind), we find 
that the steady-state accretion velocity in the layer is related to the the electron density by 

^ inene dz 

On the other hand, local angular momentum conservation requires 

Vr — i^lr ) = (10) 

ar AiTp oz 

Eliminating B^ and Vr between these two expressions leads to 

_ cnp 

Xe = ~ , (11) 

uh 2eB^nH 

where uh is the density of hydrogen nuclei, so that p/nn ~ lArUp at cosmic abundance. 
Taking B^ = —(3/4)5^5^ = 3MQ/8r, as is true for the minimal configuration ([7]), we have 

Xe ^ 8.3 X 10~^^ MZj^\~l/* , (12) 
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which is independent of Sq. Note that B • fl < in order that a;e > in eq. f[TT|) . Strictly 
steady wind-driven accretion in the extreme Hall limit ^ 1 ^ /5e) appears to be a 
degenerate case. It may be that the active layer for wind-driven accretion mus t be in the 
ambipolar (1 < A ^ Pe) regime, as in the models of IWardle fc Koenigll (Il993l ). Still, the 
electron fraction (fT2l) is interesting as a characteristic value. 

A lower bound on can be obtained by considering that the gas pressure at the base of 
the layer can be no smaller than the change in magnetic pressure across the layer. The former 
is -Pgas = ^aC^/ha ~ 2'EaCsfl, uuless the field is strong enough to significantly c ompress the 
layer (which is actually a requirement of the models of lWardle fc Koenigll Il993l ). The latter 
is APmag = (-Br + -B|)/87r, which becomes (5/8)-B^/87r in the minimum-energy magnetic 
configuration (^^. Hence 

S„>0.061M_7r;;^/'gcm-2. (13) 



This is well into the ambipolar regime, because the ion Hall parameter /3j > 1 for S < 
if the minimal field ([7]) is used for B in eq. ([8]) . 



Q.6M]_{^g cm-2 



2.3. Magnetorotational turbulence 

Accretion may be driven by radial transport of angular momentum within the disk or 
active layer rather than the torques exerted by a wind. We assume that the transport is due 
to turbulence excited by the magnetorotational instability (hereafter MRI). In a statistical 
steady state, the constancy of the angular momentum of the disk within radius r requires 

oo 

Mflr^ + 2nr^ j (p^ - B^B^/An) dz = To. (14) 

— oo 

Here Fq is the torque exerted on the inner edge of the disk; at radii 0.1 AU, it is probably 
safe to neglect this torque compared to the larger terms on the left hand side above. The two 
terms in the integrand are the relevant components of the Reynolds and Maxwell stresses 
for angular- momentum transport, with v' being the mass- weighted velocity fluctuation as- 
sociated with the turbulence, v' = v — v so that v' = 0. Linear analysis and nonlinear 
simulations indicate that the Reynolds stress i n MRI turbulence has the same sign as the 



Maxwell stress but is smaller by a factor ~ 4 (iPessah et al.ll2006l . and references therein), 
so for the purposes of the estimates below, the Reynolds stress will be neglected, and the 
Maxwell stress will be assumed to be confined to active layers each of column density 
and thickness ha on either side of the disk midplane. 
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With these approximations, the average Maxwell stress within the active layers becomes 



(15) 



Since + Bl > 2\BrB^ 
the active layers is 



a lower bound for the root-mean-square total magnetic field in 
E>3.2M!(V-"/'G, (16) 



in which we have taken the thickness of the layer to be ha ~ 0.5cs/Q as before. At this field 
strength, the magnetic pressure within the layer is comparable to the gas pressure: 



gas 



mag 



52 



7 ' AU 



;i7) 



The ratio f[T7j) is often denoted by /3 in the MRI literature but is not to be confused with the 
H all parameters f ^2.2p . The linear analysis of Wardle (1999) and the nonlinear simulations 
of ISano fc Stond (120021 ) indicate that MRI may grow in stronger fields and produce stronger 
turbulence with the Hall terms than without them, at least for the favorable sign of -B • f2. 
However, it is generally presumed that neither linear MR I nor MRI-driven turbulence can 



operate unless Pgas > Pmag (e.g., iKim &: Ostrikerl l2000l ). so eq. (fTTl) suggests that > 



10M_7 r~l g cm~2 if MRI dominates the accretion process. The ratio of the minimum fields 
( IT6ll and (JTj) is ~ ^/rjh^, because the turbulent stress is exerted over an area proportional 
to the thickness of the active layer, whereas the wind stress is exerted on the (much larger) 
horizontal surfaces of these layers. 

Inserting the lower limit flTBl) into the expressions for the Hall parameters yields 



A>0.61Sr^M 



1/2, 



, so the ions may be marginally magnetized. For /5j > 1, the slippage 



between the neutrals and the field lines is further enhan ced by ambipolar diffusion, which 
scales oc B^ (IWardle fc Koenigllll993l : IWardlelll999l . 120071 1. This is somewhat complicated to 
discuss if the abundance of negatively charged grains is comparable to that of free electrons. 
In view of the discussion of magnetic and gas pressures above, however, we expect that (3^ 
will remain less than unity — though not by much — at the radii and accretion rates of interest 
to us, so that the Hall diffusivity is still marginally dominant, riH = (? /^TicfH ~ Bc/ATreUe. 
(A more accurate treatment would likely lead to larger lower bounds on Xe-) A dimensionless 
measure of the coupling between the neutrals and the magnetic field within the active layer 
is the magnet ic "Reynolds number" based on this Hall diffusivity, ReM,H = Csha/rju] hke 
Wardld (119991 ) . we presume that ReM,H ^ 1 is necessary for sustained turbulence. In that 
case. He > cBQ/27iecj., and therefore 



X, 



> 2.6 X 10""My7^ S^V 



-9/8 



cm 
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As noted above, the ratio of Ohmic to Hall conductivities is \(yo/<^H\ ~ I3e when /?e S> 1 
and (3i < 1. At the lower bound ^ for the field strength, (3^ ~ 300My/S^Vi{f . It follows 
that the magnetic Reynolds number f[T^ based on the Ohmic diffusivity should satisfy 
Rcm ^ 100 in order to support MRI-driven accretion at rates > lO'^MQyr"^. 



Sano fc Stond (120021 ) and [Turner et al.l (120071 ) found that a criterion based on the El 



sasser number, A = v\^/ (r/o^) = 1; accurately traces the edge of the MRI turbulence in their 
simulations. (But the former authors referred to this quantity as "magnetic Reynolds num- 
ber".) The Elsasser number depends explicitly upon magnetic field as well as Xg. Using A to 
define wo uld impede comp a rison -w ith most previous studies of the ionization balance — 
in particular, lllgner &: NelsonI (l2006al ) — which have adopted thresholds in c^/rjofi. In fact, 
however, linear stability depends upon at least two dimensionless parameters: (Rcm,-^), 
or equivalently (c^/V^^, A) since Re^/A = c^/Vl^ The criterion Re^ > 100 for the active 
layer is more conservative (i.e., allows a larger active layer) when c^/V|^ > 100, which is the 
case in almost all published MRI simulations because it is desired that the fastest-growing 
vertical wavelength b e smaller than the gas scale height Cs/Q. For example cj./'yV^z > 400 in 
Sano fc Stone J2002h. where 7 == 5/3 is the adopted adiabatic index. Moreover, this criterion 
is consistent with lTurner et al.l (120071 ) 's result (see Fig. 6 of their paper). A strong — but not 



too strong — vertical field may permit linear growth at Rcm ~ 1- iGammie fc Balbus I (119941 ) 
found instability up to Vaz/cs < 1.5 in their "quasi-global" analysis of a vertically stratified 
shearing box, although their analysis was limited to ideal MHD. The nonlinear development 
of MRI with equipartition-strength background fields has not been well explored, but it is 
conceivable that radial advection of flux somehow prefers such states. 



2.4. Summary of the constraints on field strength and ionization fraction 

For M ~ lO^^M0yr^^, wind-driven and MRI-driven accretion require similar ionization 
fractions, Xe ~ lO^^^-lO^^'^cm"^, as shown by equations ( fT2l) and ( fTSj) . Winds may allow 
weaker fields than turbulence for the same accretion rate [eqs. (I7j) vs. ( |T6i) ] because of a 
geometric advantage: wind torques are exerted on the disk surface, which has a larger area 



than the disk thickness by ~ r/h. However, winds have their own theoretical difficulties (jShu 



I991I : lOgilvie fc Livid l200ll ). and the field may have to be much stronger than the minimal 
value (^^. The constraint that the field not exceed equipartition allows winds to couple to a 
smaller surface density than turbulence [eqs. (|T3ll & (ITTl)]. 

We emphasize that the constraints on field strength, active surface density, and ioniza- 
tion fraction discussed in §2.31 are based on accretion rates in the nonlinear regime rather 
than the minimal conditions for linear MRI instability, which are more easily satisfied. 
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3. Conductivity Calculation: Methods 

In this section we describe our model for the disk conductivity, including sources of 
ionization, chemical reaction rates, and grain interactions. Readers interested only in the 
results may wish to skip ahead to ^ 

Our criterion for the active layer is Rcm > 100, where 

ReM = ^ , (19) 

is the magnetic Reynolds number, and 

rio = 230X^/23,-1^^2 g-i ^20) 



is the magnetic diffusivity based on the Ohmic conductivity ao ~ e rie/ me{crv)en (iBlaes fc Balbus 



19941 ) ■ The discussion of ^indicates that the use of the single dimensionless parameter (|T9l) 
is an oversimplification — the Hall parameters ([H]), which involve the field strength, are also 
important — but the criterion Rcm > 100 appears to be necessary if not sufficient near 1 AU, 
at least if MRl turbulence dominates the transport of angular momentum. 



3.1. Ionization Model 



A number of nonthermal processes may contribute to an excess of free electrons over 
the abundance in LTE. We consider X-rays from the protostar and cosmic-ray ionization. 

T Tauri X-rays can be very energetic, wi th luminosities Lx ~ 10^^ — lO^^erg s~^ and pho- 



ton energies ranging from about 1 to 5 keV (jCasanova et al.lll995l : ICarkner et al.lll996l ). We 
model the X-ray source by two bremsstrahlung-emitting coronal rings at r adii 10-Rq from the 
rotation axis and a similar distance above and be l ow the disk midplan e (IKrolik fc Kallman 
19831 : Iciassgold et"aDll997l : llgea fc Glassgoldll 19991 : IPromang et aDl2002h . The X-ray photons 
are attenuated by photoionization (absorption) and Compton scattering. Scattering reflects 
some of the X-rays photons from the disk, therefore reducing the total ionization rate; but 
it also allows photons to penetrate deeper into the disk by deflecting oblique rays towards 
normal incidence. The photoionization cross section for keV X-ray photons is ~ 10~22cm2, 
decreases roughly as E~^^^, and falls below the Thomson scattering cross section at -Ephot ^6- 
7 ke V. We take the X-ray ion ization rate as a function of column density normal to the disk 
from llgea fc Glassgoldl (119991 ). who performed Monte-Carlo radiative transfer calculation in- 
cluding scattering in the MMSN model, assuming that metals are depleted onto grains and 
segregated from the gas. For Tx = 3keV, their results imply an ionization rate per hydrogen 
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molecul^ can be well fit by 

Cf ( R 



lAU 



2.2 



(21) 



where Lx,29 = Lx/lO^^erg s ^, A^hi,2 is the column density of hydrogen nucleus vertically 
above and below the point of interest, R is the cylindrical radius to the central protostar, 
Ci = 6.0 X 10-i2g-i^ ^ X l02icm-2, a = 0.4, C2 = 1-0 x 10-^5s-\ N2 = 7.0 x lO^^cm-^, 
(3 = 0.65. The first exponential represents attenuation of X-ray photons by absorption, while 
the second exponential incorporates a contribution from scattering. The scaling with radius 
is slightly steeper than inverse square but less steep than it would be without scattering. 
For Tx = 5keV, we fit their results with (i = 4.0 x lO^^^g-i^ jy^ ^ 3 9 x lO^^cm'^, a = 0.5, 
C2 = 2.0 X 10"i^s-\ N2 = 1.0 X lO^^cm-^, p = 0.7. We have also compared this fitted 
ionization rate to a direct calculati o n of X -ray ionization rate without scattering based on 
equations (2)-(4) of iFromang et al.l (120021 ) (see Appendix A for supplemental information). 
The latter gives a slightly larger ionization rate at 1 AU. However, at larger radii (e.g. 
10 Au), scattering increases the ionization rate at columns S > lgcm~^. Throughout this 
paper, unless stated otherwise, we take kTx = 3keV and Lx = 5 x lO^^ergs"^. The 



ionization rate depends rather weakly on the X-ray temperature in llgea fc Glassgoldl (Il999l )'s 
results [see the i r Fig. 3], and much more sensitively with the pure-absorption formalism of 
Fromang et al. ( 2OO2I ). We do not entirely understand these differences, but w e base the 
calculations described below on our fits to the results of llgea fc Glassgoldl (119991 ). 



Interstellar cosmic rays (CR), unless shielded by the stellar wind, provide ionization 



Aeff 
S.CR 



Coexp(— S/96g cm 



(22) 



per hydrogen atom with stopping grammage Eg ~ 96g cm ^ (jUmebayashi fc Nakandll98ll ). 



Recent ob servations of the co smic-ray fiux towards the diffuse cloud ( Persei suggest (0 ~ 
10^^^ s"^ (iMcCall et al.ll2003l ). an order of magnitude larger than older values. We find that 
cosmic rays have little effect on for r < IOau unless Co ^ lO^^^s"^, so we adopt Co = 
as our standard value but also consider values up to Co = 10^ 



We neglect photo-ionization by ultraviolet photons because of their small penetration 
depth. We also neg lect the thermal ionization of the alkali metals Na"*" and (e.g., 
Fromang et al.ll2002l ). which becomes important only above ~ lOOOK, and ionization by 



radionuclides, which appears to be negligible. 



^We adopt this definition of ionization r ate throughout th i s pap er, which is 2 times larger than the 
ionization rate per hydrogen atom. Note that Igea fc Glassgold ( 1999[) adopt the latter definition. 
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The total effective ionization rate now is ( = (x + Ccr- Hydrogen and helium are the 
main targets, so we include only the four ionization reactions listed Table [1] in our chemical 
network. Ionization of atomic hydrogen, neglected by IN06a, is included because in the 
absence of grains and therefore of H2 formation, chemical equilibria involve H rather than 
H2. The third reaction ensures that the ionization rate is almost independent of the division 
between atomic and molecular phases (see §3.4p . 

X-ray ionization is much stronger near the star, but the X-ray flux falls as and drops 
rapidly toward the disk midplane. As a result. X-rays ionization dominates at modest radii 
(r < 5 Au) and near the surface. Cosmic ray ionization is almost negligible in the inner parts 
of the disk but dominates in the outer disk and toward the midplane. 



3.2. Chemical Reactions 



We consider two alternative chemica! 



introduced by lOppenheimer fc Dalgarnd ( 



reacti on networks. The first is a simplified network 



Gammie 



9741) for dense molecular c l ouds but also widely 



19961 : iGlassgold et al.lll997l : iFromang et al.ll2002 



applied to protoste llar disks (e.g 
Turner et al.l 120071 ). This network is basically a two-element kinetic model involving five 



species: a molecular species m; a neutral atomic gas-phase metal M; their ionized counter- 
parts m"*" and M"*", and free electrons e~. There are four reactions, as listed in Table 1 of 
IN06a, whose reaction rates we adopt. 

Following IN06a, we also consider a much extended kinetic model involving nine elements 
(H, He, C, O, N, S, Si, Mg, and Fe) and 174 species as given by Table A.l of IN06a. We 
extract 2109 reaction s involving these species from the latest version of UMIST database 
( IWoodall et al.l 120071 ). In our selection, we neglect photo reactions because of the poor 
penetration of UV into the disk. We also neglect "collider" reactions, which have very low 
rate coefficients. We further exclude reactions involving CR protons and CR photons in the 
UMIST database, since we have already included the dominant reactions shown in Table [TJ 
There are in total 2113 gas-phase chemical reactions including 4 ionization reactions in Table 



1. 


H2 - 


H++e- 


0.97 


^eff 


2. 


H2 - 


H+ + H + e- 


0.03 


^eff 


3. 


H - 


H+ + e- 


0.50 


^eff 


4. 


He - 


He+ + e- 


0.84 


^eff 



Table 1: Ionization reactions. Final column is the ionization rate for the reaction shown. 
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[U We have 147 more reactions than IN06a, who used an older version of the database q 

The rate coefficients in the database are functions of temperature. Where the local disk 
temperature lies outside the stated range of validity, we adopt the same procedure as IN06a: 
we replace Tdisk with the upper or lower bound of the vahd range, as appropriate, before 
computing the rate. 



3.3. Reactions with Dust Grains 

Dust grains are efficient absorbers of free electrons. Grains also interact with other 
neutral and ionized species mainly by adsorption, desorption and charge exchange. Further, 
adsorbed species hop on grain surfaces, and reactions can take place on grain surfaces. There 
is currently no standard database for grain reactions. We follow the prescriptions of IN06a 
involving "mantle chemistry" and "grain chemistry" . We assume that the maximum charge 
of a grain particle is 2, which is appropriate for very small Xe- 

Mantle species are adsorbed counterparts of gas-phase neutral species. Reactions that 
involve mantle species belong to mantle chemistry, as given by Table 3 of IN06a. A mantle 
species X[m] is formed by collisions between a neutral species X or its ionized counterpart 
X"*" with a grain particle, and is destroyed by desorption. The mantle species defined here 
are independent of grain charge. 

Other grain-related reactions belong to grain chemistry, as given by Table 4 of IN06a. 
These include charge-exchange reactions between ionized species and negatively charged 
grains, absorption of free electrons by grains, and grain charge-exchange reactions. Several 
ionized species do not have any neutral counterpart (e.g. Hg , HaO"*"), as needed in the 
charge-exchange reactions. Following IN06a, we assume the products of these reactions to 
be the same as for dissociative reactions in the gas phase, e.g. Hj]" -|- gr~ — > 3H -|- gr. If there 
are multiple final states, we adopt the gas-phase branching ratio. 

Below we describe our adopted grain properties and outline the calculation of rate 
coefficients for grain reactions. 



^Note that in the new version of UMIST database, several species names are changed. 
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3.3.1. Assumptions about Grains 

All grains are taken to be spherical with density pd = 3gcm"'^. Ideally, one should 
consider a size distribu tion of grains . The conventional choice for interstellar grains is the 
MRN size distribution jMathis et al.lll977h 



N{a)da oc a da, < a < amax , (23) 

where a is grain radius, and the usual cutoffs are amin ~ 0.005/im and amax ~ 0.25/im. In 



many previous works on disk ionization (e.g., Illgner fc Nelson! l2006al : ISalmeron fc Wardle 



20081 ). grains are chosen to have a single size. We allow up to two grain sizes, oi and 02, 
with different abundances. As standard values, we take ai = 0.01/im and 02 = O.l/im. The 
mass ratio of the two grain populations should be /1//2 = (01/02)"°'^ to crudely represent 
the MRN distribution (1231) . with / = /i + /2 = 0.01 as the standard mass fraction. The 
grains are taken to be uniformly mixed, though we vary both the sizes and the abundances 
of the two populations to imitate grain growth and precipitation out of the surface layers. 

Each grain population has its own mantle and grain chemistry. Collisions between the 
two populations are allowed and cause charge exchange in the same way as collisions among 
members of the same size population. We assign equal branching ratios to the transfer of 
one or two electrons. In fact, we find that collisions between grains belonging to different 
size populations are unimportant. We will show that the electron abundance, and therefore 
the size of the active zone, is dominated by the smallest grains. 



3.3.2. Grain reaction rates 

1. Collisions between ions/electrons and grains. These reactions correspond to reac- 
tions 7, 8 in Table 3 and reactions 1-6 in Table 4 of IN06a. They recombine free electrons and 
add to the Ohmic resistivity. The rate coefficients involve collision rates and the probability 
that an electron stays on a grain after colliding with it (sticking coefficient). 



The collision rate is estimated following equations (3.1) and (3.3)-(3.5) of lDraine fc Sutin 



( 119871 ). it is modulated by Coulomb interactions and induced polarization. The sticking co- 
efficient S is estimated as follows. For ion-grain collisions 5* is insensitive to temperature, 
so we take Sx+ = 1- The sticking coefficient Sp, for elec t ron-io n collisions is more subtle. 



We adopt equations (B5), (B6), and (B13) of iNishi et al.l (Il99ll ). Electrons undergo elastic 



and inelastic scatterings (assumed to be isotropic) with grains until absorbed. The elec- 
tron sticking coefficient depends sensitively on temperature (IN06a). For the purpose of 
calculating Se, we take the grains to be made of graphite, with atomic weight 12, Debye 
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temperature 420K, and electron binding energy = 1 eV, though 5*6 is insensitive to the 
latter parameter. 

2. Collisions between neutral gas-phase particles and grains. These correspond to 
reactions 1-5 of Table 3 in IN06a. These reactions can potentially deplete gas-phase species by 
adsorbing them onto grain surfaces, if the inverse of this process (desorption) is slow. Similar 
to the previous case, the rate coefficients are expressed by the product of collision rate and 
sticking coefficient. The collision rate is just geometric. For incident particle energy Ei, the 
probability of being adsorbe d is approximately = exp(— e^/2), where e = Ei/ y/WKEl 



( iHoUenbach fc Salpeterlll970l ); D and Ai?s denote the dissociation energy and the amount 
of energy transferred to the grain particle as lattice vibrations, respectively. For each neutral 
gas-phase particle, D is approximated by its binding energy E£, as given in Table A 2 of 
IN06a, and AE's is approximated by 2.0 x 10^^ eV. The sticking coefficient is then obtained 
by integrating over a thermal energy distribution of Ei. 



3. Collisions between grains. These are reactions 7-10 in Table 4 of IN06a. They 
re distribute charge among grain particles. We calculate the collision rates from equation (3) 



of lUmebayashi fc Nakand (jl990l ). Note that we take pd = 3 gem ^. Since the abundance 



and thermal velocities of grains are low, these reactions occur very slowly. 

4. Desorption processes. These are described by reaction 6 in Table 3 of IN06a. 
Mantle species migrat e among surface sites at a characteristic thermal hopping frequency 



(jHasegawa et al.lll992l ) 



V TT^m 

where Us ~ lO^^cm"^ is the surface number density of binding sites, and m is the mass of 
the adsorbed species. The desorption rate of mantle species via thermal hopping is then 

k = uoexp{-ED/kTd) , (25) 

where is the temperature of the dust, taken to be the same as the gas temperature. 

Note that both adsorption and desorption rates have an exponential factor. For adsorp- 
tion rate, the exponential factor is ~ exp [— (T/Tq)^], where Tq depends on Ej:, and grain size. 
For desorption rate, the exponential factor is exp {—Ed/T). At relatively high temperature 
(e.g. T > lOOK for metals), the desorption rate is much higher than adsorption rate (since uq 
is quite large), hence almost all species are in the gas phase. As the temperature decreases. 
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the mantle abundance increases super-exponentially. This is particularly pronounced for 
metals because they have large binding energy E£,. We give an example to demonstrate this 
effect, which will aid in interpreting the results in §4.21 

The binding energy of magnesium is En = 5300 K. Balancing adsorption and desorption 
processes, we obtain 

nwg _ 1-34 X 10^5 exp (-5300K/T) exp [(T/496K)2] 
riMgH ~ (ngr/lcm-3) x (r/300K)i/2(a/i^ni) ' ^ ' 

For Ugj- = O.lcm^^, a = O.l/im, and T = 300K, the ratio (12^ is about 4.1 x 10^. However, at 
lOOK, the ratio becomes as small as 2.3 x 10~^. The critical temperature T ^ 150K, below 
which most metals are adsorbed onto grains. 



3.4. H2 Formation on Grain Surfaces 

In the interstellar medium, grains catalyze many reactions. Adsorbed atoms and molecules 
hop among sites until they collide, and the heat of reaction is absorbed by grain lattices. 
Most importantly, H2 is very efficiently formed on grains in dense molecular clouds, where 
T < 20K. Whether H2 formatio n on grains is still i mportant a. t the higher temperatures 
of protostellar disks is unknown. ICazaux fc Tielend (120021 . |200J) proposed a model incor- 
porating both physisorption and chemisorption of H atoms on grains. Physisorption is the 
adsorption process mentioned above, whereas chemisorption involves much stronger (~ 1 eV) 
chemical bonds. Radical species with unpaired electrons are likely to be chemisorbed. At low 
temperatures, H2 forms mainly by interactions between a physisorbed and a chemisorbed H 
atom; at high temperature, two chemisorbed H atoms are involved. The works cited above 
found th at H2 formation can be efficient up to about 500K, with efficiency up to about 0.2. 
However, Cazaux fc Tielens (2002, 2004) overestimated their H2 formation efficiencyjfl 



Our chemical evolution model shows that if H2 formation is not considered, almost all 
hydrogen will ultimately be converted into atomic form, an unlikely state for a protostellar 
disk. Therefore, we include H2 formation on grains, but with low efficiency. In view of the 
physical uncertainties, we model this process by assuming a uniform probability t] for a pair 
of adsorbed (physisorbed) hydrogen atoms to form an H2 molecule. Then the reaction rate 



■^Equations (1) and (2) of Cazaux fc TielensI (2004) omit a factor ^{E — Bij)/E, thereby overestimating 
the transition rate from physisorbed sites to chemisorbed sites and overestimating the inverse transitions. 
Hence the surface number density of chemisorbed H should be lower than they claim. Wc have confirmed 
this with the author (Cazaux, private communication). Note that there are other typos in the formula. 
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for H + H ^ H2 is 

R = ^nnVHridadiiSH , (27) 

where uh, Ud are the number densities of H atoms and grains, vh is the atomic thermal 
velocity, ad = t^o? is the geometric cross section of grains, and Sh is the sticking coefficient. 
For all 7] > 10~^, hydrogen remains molecular. We take rj = 10~^ as our standard value. 



3.5. Initial Conditions 

We normalize the abundance (proportional to number density) of hydrogen atoms to 
unity, and all other elements and grains proportionately, with the relative elemental abun- 
dances given in Table 6 of IN06a. Elements in the gas phase other than hydrogen and 
helium are depleted compared with solar abundance. We vary the metal (Mg, Fe) and grain 
abundances. Note that we use Xi to denote elemental abundances, and x[X] to denote the 
abundance of species X. All grains are initially neutral, and all elements in their atomic form 
except hydrogen. 



3.6. Numerical Method 

The chemical evolution equations are a set of first order ordinary differential equations 
(ODEs). Due to a broad ranges of reaction rates, the se equations are very stiff. We use 



the stiff . c subroutine described in iPress et al.l (119921 ) as the main integrator. This uses 



the Kaps-Rentrop algorithm, a 4th order implicit method. The evolution equations conserve 
charge and elemental abundances, but numerical errors violate the conservation laws. In our 
program, we enforce conservation after each step by adjusting the elemental abundances and 
Xe- The integrator efficiently evolves the system for 10^"'' years. 

We calculate the electron abundance by evolving the system until chemical equilibrium is 
reached. Absent irreversible reactions, chemical equilibrium would be unique and would obey 
detailed balance. The inverses of about 70% of the reactions in our network are neglected, 
however, because these inverses are too slow to be effective. For example, as three-body 
reactions are very slow at the low densities of interest, reactions with more than 2 products do 
not have any inverse in the reaction chain. Because there are so many irreversible reactions, 
it is conceivable that the equilibrium state is not unique. We have tested a few different 
initial conditions, and the abundances of the main species converge to the same values when 
integrated long enough. Some radiative association reactions can be so slow that equilibrium 
is approached only after very long time. Since protostellar disk lifetimes are believed to be 
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at most a few million years, we evolve the equations for 10^ years and accept the final Xe as 
the equilibrium value (except in §4.11 where 10^ year is used)B 

We determine S^, the column density of the active layer, as follows. We assume that 
Rcm increases monotonically upward. If the midplane is active, Reuiz = 0) > 100, then 
the whole column is active. Otherwise, we start at a large disk height where Rcm > 100 and 
search by bisection to find the height at which Rcm = 100. 

Our model parameters are summarized in Table [2] together with their standard values 
and ranges. For a single population of grains, a = O.l/im and / = 0.01 are the default 
values. Not listed in the Table are some fixed parameters mentioned in the text, such as the 
power-law indices of T(r) and S(r). 



4. Conductivity Calculation: Results 

In the following subsections, we start from grain- free models ( §4.11) . and gradually in- 
crease model complexity by adding one ( §4.21) and then two ( §4.31) populations of dust grains. 
In each subsection, we first discuss the chemistry, evaluating the free-electron abundance as 
a function of density, temperature, ionization rate, metal abundance, and grain properties (if 
applicable). Then we apply these results to the MMSN disk model to calculate the thickness 
of the active zone and its dependence on model parameters such as X-ray luminosity, metal 
abundance and grain properties. In both steps, we compare the simple chemical network 
with the complex network. 



4.1. Models without Dust 

As a first step, we study the free electron abundance in pure gas-phase chemical reaction 
networks. In general, this abundance Xe is a function of density (p), temperature (T), effective 
ionization rate (C^^), and metal abundance (xMg and xpe)- In the simple model, the metal 
abundance is just the abundance of Mg. In the complex model, we use the combination 

Xm — ^^Mg "I" ^^Fe with XMg = 4xFe- 

Our network runs very efficiently for pure gas-phase chemical reactions. It takes less 
than 1000 years for the simple network to reach chemical equilibrium. For the complex 
network, we find the electron abundance is still subject to very slow variation after evolving 



''Even 10^ yr is long compared to the time spent by most gas elements at constant p, T, and C- 
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for 10^ years. Although 10^ year is somewhat longer than the lifetime of the accretion phase 
of the protostellar disk, this is balanced by our artificial choice of purely atomic initial species. 
Therefore we choose 10'' year as the standard evolution time in this subsection only. 



4-1.1. Chemistry 

In Figure [1] we plot the free electron abundance as a function of gas density and tem- 
perature for both simple and complex chemical models. We also vary the effective ionization 
rate and metal abundance in each plot. 

Note that gas densities within the disk span several orders of magnitude. Since the 
(unshielded) ionization rate per unit volume is proportional to gas density, and the (two- 
body) recombination rate to density squared, one expects Xe oc p~^/^. In Fig. [1^, the slope 
is slightly flatter, Xg oc p~°'^. By the same token, since ionization is almost the only source 
of free electrons, one expects Xe oc {C^Y^'^. In both panels of Fig. [H one sees that as 
varies by a factor of 100, Xe varies by ~ 7-10. 

One sees also that Xe is insensitive to temperature. This is because ionization reactions 
are independent of T, while recombination reactions are exothermic, with rates roughly 
proportional to T~^/^. Most other reactions scale between T~^/^ and T^/^. In protostellar 
disks, where T ~ 100 — 1000 K, T^^/^ varies much less than density and ionization rate and 
therefore affects Xe only slightly. 

In these grain-free models, the electron abundance is very sensitive to metal abundance. 
Metal atoms donate their electrons to ions in charge-exchange reactions (mainly with 
and H^)] the resulting metallic ions recombine with free electrons only by radiative reac- 
tions, which are much slower than the dissociative reactions available to molecular ions, e.g. 
H+ + e- > 3H. 

In most circumstances, the electron abundance obtained from the complex network is 
greater than that obtained from the simple network by roughly a factor of two. Here we 
differ from IN06a, who found the opposite tendency. This i s due in part to our use of a 
later version of the UMIST database for the complex network( Woodall et al. 2007 1^. There 



are notable changes of reaction rate in a number of reactions in the new database, and in 



^Also, in the absence of dust, the complex network graduaUy converts molecular to atomic hydrogen (see 
section [3.4[) . whereas the simple network subsumes hydrogen into the molecular species m and m+. Our 
addition of H ionization (see section [3. 1[) increases the ionization rate. If we remove H ionization, electron 
abundances for two networks are similar. 
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chemical equilibrium the abundances of several species are also quite different. Recently, 
Vasyunin et al.l (120081 ) performed a sensitivity analysis of the UMIST06 database. It was 



found that typical uncertainties of molecular abundances do not exceed a factor of 3-4. 
The differences in Xg between our results and those obtained by IN06a with the UMIST99 
database are comparable to these uncertainties. 



4.1.2. Application to the Disk Model 

Now temperature and density are fully fixed by the disk model, and our free parameters 
are ionization rate and metal abundance. The ionization rate is characterized by 3 parame- 
ters, namely. X-ray luminosity (Lx), X-ray temperature (Tx) and cosmic-ray ionization rate 
Co- Since we already know how metal abundances and ionization rate can affect Xg, and we 
also know how ionization parameters affect ionization rate in the disk (see §3.10 . our main 
goal now is to check how these parameters affect the column density of the active layer, 
(Note that is defined as the active column on one side of the disk). In Figure [2] we plot 
Sa versus radius. We also compare simple and complex chemical networks in Fig. [2l 

Figure [2] shows that in the absence of dust grains, the requirement > lOgcm"^ (§2) 
is easily realized. Beyond 2au, the whole disk becomes active once cosmic-ray ionization 
with > 10^^''s~^ is turned on. Without cosmic-rays. X-rays can also render the whole disk 
active beyond 6 AU. In the absence of dust grains, the metal abundance strongly influences 
the thickness of the active layer. At 1 AU, reduction of xu from 1.25 x 10"*^ to 1.25 x 10^^^ 
causes to drop by a factor 10 for the complex network (20 for the simple network). Also, 
we see that the complex model produces a slightly thicker active zone, consistent with the 
result in §4.1.1[ 



4.2. Single-sized Grains 

In this subsection, we add a single population of dust grain particles to the chemical 
reaction network. The standard grain size is a = O.l/xm, and the standard mass fraction is 
/ = 0.01. We evolve the reaction network for 10^ years. 



4.2.1. Chemistry 



In Figure [3] we plot the free electron abundance as a function of gas density for simple 
and complex chemical reaction networks respectively. As before, we find that Xe scales with 
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Fig. 1. — Electron abundance Xe in grain-free models, (a) Xe versus gas density p at constant 
T — 280K. (b) Xe versus T at constant p — 10~^^gcm~^ . Dash lines: simple network. 
Solid lines: complex network. Bold lines are for standard parameters: C,^^ — 10~^^s~^, 
xyi — 1-25 X 10~^. Each of the other curves differs from the standard in one parameter, as 
shown. 
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a) Simple chemical network 




1 10 50 

Radius (AU) 



b) Complex chemical network 




Radius (AU) 



Fig. 2. — The column density of the active layer defined by Rcm > 100 (§2) versus radius 
for simple (a) and complex (b) chemical reaction networks without grain particles. The 
upper bold line indicates the half the total surface density of the disk. The middle bold line 
corresponds to our standard parameters: Lx = 0.5 x 10^°ergs~^, Tx = 3keV, Co = Os~^, 
xm = 1-25 X 10~®. Other curves differ from the standard curve in one parameter, as marked. 
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a) Simple chemical network 
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b) Complex chemical network 
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Fig. 3. — Electron abundance versus gas density at constant temperature, T — 280K for 
simple (a) and complex (b) reaction networks with a single population of grains. The bold 
lines are for the standard values C,'^^ — 10""^''' s~^, xm — 1-25 x 10~^, a = O.l/xm, / = 0.01, 
T] — 10^^. Other curves differ in one of these parameters, as marked. Bold dashed lines are 
grain- free (/ = 0). Light dashed line in (a) is for reduced metals, as shown. Note the tiny 
difference. Panel (b) is even less sensitive to xm- 
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Fig. 4. — Electron abundance versus grain size with a single population of grains. Solid 
lines: complex network; dashed hues: simple network. Plots with bold lines have fixed total 
grain area f/a =const, while plots with thin fines keep f/a^ constant. Panel (a): Bold 
lines: f/a = OM/im'^; thin lines: f/a^ = 0.01//m~^. (b) Bold lines: f/a = lO^^iim"^; 
thin lines: f /a^ = 10~^/xm"^. Lines in the upper, middle and bottom groups correspond 
to p = 10-l^ 10-^^ and IQ-^gcm-^ respectively. In aU plots, T = 280K, = 10-^^s-\ 
xm — 1-25 X 10~^, T] — 10~^. Note that in panel (b) the grain size is extended to 10//m 
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density and ionization rate, and is not very sensitive to temperature. There are also some 
notable differences. Firstly, comparing the bold solid and dashed lines, one sees that with 
dust grains, decreases faster as density increases, suggesting that dust is more effective 
in suppressing ionization in denser environments. Secondly, comparing Fig. [3] with Fig. [T^, 
one sees that changing ionization rate by a factor of 100 causes Xe to change by less than a 
factor of 10 in the absence of dust, but by a factor of 25 in its presence. 

A striking difference from the grain-free case is that gas-phase metal atoms are no longer 
important. This is consistent with what IN06a have found, but our explanation is somewhat 
different from theirs. At 280K, metal atoms are not swept up by grains (see §3.3.21) at all. 
Rather, the role that grains play is to promote recombination. We find that when we add 
grains, the number density of metallic ions is greatly reduced. 

With our fiducial parameters, we see that the complex network produces more free 
electrons, by a factor of < 2. However, as we gradually suppress grains, the simple model 
recovers the grain-free result more rapidly than the complex model does. In certain ranges, 
the simple network can produce larger x^ than the complex one. For the simple network, 
reduction of grains by lO'' (/ = 10^^) is almost enough to recover grain-free result in low 
density regions, but for complex network, depletion of nearly 10® (/ = 10^*^) is required. 

In Figure [31 we see that electron abundance increases as grains are suppressed, and 
at fixed grain mass fraction, Xe increases with grain size. This is as expected, but grains 
actually take effect in a complicated way, especially in the complex network, as suggested 
by the fact that the curves in Fig. [3] are not simply vertical translates of one another. It is 
natural to ask what combination of a and / best controls Xe- To investigate this, we plot in 
Figure H] Xe versus grain size a for a single population of grains at various densities. Plots 
with bold lines have fixed total surface area No? (x f /a =const, while for thin lines we have 
fixed Na oc //a^ =const. If total surface area completely determines Xg, we would expect 
that bold lines to be fiat. In Fig. H^, the bold lines for both complex (solid) and simple 
(dashed) network are approximately flat, and Xe slightly decreases as a decreases. This 
means that total grain surface area is a good approximation to the controlling parameter, 
though small grains are slightly more effective than big ones. In Fig. HJo, as grains are 
substantially depleted, we see more complicated dependencies, but with fixed total surface 
area, the trend is still that Xe decreases as a decreases. The fact that small grains are more 
efficient in reducing Xe at fixed total surface area suggests that the controlling parameter may 
lie somewhere between the total surface area Na"^ oc //a and the combination Na oc f /a^. 

The results above can be qualitatively understood. The cross sections of all grain 
reactions scale with grain surface area, but they are modulated by a grain-polarization term 
that enhances the cross section of smaller grains. For the reaction e~+gr gr~, the collision 
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rate per grain is proportional to (jPraine fc SutinI 119871 ) 



^'^[df) =^^0.30[—^ (28) 

and for e~ + gr+ — >■ gr, the reaction rate is roughly proportional to 

S,.^^^.^,om,?^'J^ (29) 

vra^ Zakl I a 

Although smaller grains also have smaller sticking coefficients for electrons, we find through 
numerical experiments that over a wide temperature range, the electron sticking coefficient 
for an a = 1.0/im neutral grain is greater than that for an a = 0.01/im neutral grain by a 
factor only ~ 2. Therefore, the electron abundance is roughly controlled by the total surface 
area of the dust grains, but smaller grains lead to smaller Xg at fixed total surface area. 

Our simple analysis here considers just the direct electron absorption by dust grains, 
which roughly applies to the simple network. For the complex network, however, much more 
complicated interactions between dust and all the species, as well as the complex reactions 
in the gas phase make it less predictive. Our numerical experiments above shows that in 
certain regimes, total surface area controls Xg, but in other regimes, na oc f/a^ may be a 
better approximation. 

One more comment on the two reaction networks. As mentioned above, Xg converges 
to its grain-free value more rapidly as / — > in the simple than in the complex network. 
This fact is seen more clearly by comparing the left and right panels of Fig. HI Before 
grain depletion, x^ in the complex network is typically slightly larger than that in the simple 
network (left panel). As we deplete grains by a factor of 100 (right panel), the simple network 
typically yields a larger Xe than that of the complex network, by up to a factor of 10. 



Tra^ 



4-2.2. Application to the Disk 

Figure O plots against radius in the same way as in Fig. [21 and compares different 
model parameters, especially grain size (a) and grain mass fraction (/). The column density 
of the active zone is dramatically reduced when we add sub-micron grains. Both networks 
predict < lgcm~^ at 1 AU. As in Fig. [3l is more sensitive to ionization than in 
the grain-free case. X-ray ionization alone cannot ionize the whole disk with our fiducial 
parameters; but cosmic rays, if unshielded, can ionize the entire column of the outer disk 



beyond ~ 10 AU. These features are also in agreement with ISano et al.l (120001 ). 



One sees that has a relatively simple dependence on the grain properties in Fig. [5^ 
(i.e. simple network), but the dependence in Fig. [5]d is more complicated. With fiducial 
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a) Simple chemical network 




-1 a=O.OlM,m 
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b) Complex chemical network 




1 10 50 



Radius (AU) 



Fig. 5. — Like Fig. [21 but for a single population of grains with standard values a = 0.1/im, 
/ = 0.01, and rj = 10"'^. All other parameters at standard values (Table [2]) except as 
indicated. For comparison, we plot the grain-free case in bold dash-dotted lines. 
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parameters, Sa is larger in the complex network, but with reduced dust, the simple network 
can produce a thicker active layer. Our previous conclusion that total surface area of dust 
grains roughly controls the electron abundance can be checked in the simple network, where 
we see that the a = Ifim (/ = 0.01) line is below the / = 10~^ (a = 0.1/im) line. In the 
complex network, we can see that when r < 0.8 AU, the a = 1/im line is above the / = 10""^ 
line. This corresponds to the situation of Fig. |3)d at relatively low densities. Again, the 
effect of grains on the complex network does not reduce to a single controlling parameter. 

In Fig. [5^ we see that at about 8 AU, the two solid hues for / = 10~^ and / = 10~^ 
drop sharply. This corresponds to the sudden depletion of metals onto grains at T ~ lOOK, 
as discussed in §3.3.2[ With less dust (smaller /), the transition temperature decreases 
[see eq. fl2Bl)]. and the transition radius increases. Outside the transition radius, metals are 
effectively depleted onto the grains, unless the grain number density is extremely small. In 
the latter case, our grain adsorption model may fail since the surface of all grains can be fully 
occupied by adsorbed species. Similar transitions in Fig. [5)d (e.g at / = 10~^ and / = 10^*^) 
are not well explained by the analysis in §4.2.11 

For the fiducial parameters, metals are not depleted. Under this assumption, we see 
from Fig. [5]that for to be ~ lOgcm^'^ at 1 AU, a = O.lfim dust grains must be suppressed 
by a factor 10^^ (/ = 10~^) for the simple network, and by 10~^ (/ = 10~^) for the complex 
network. Suppression by a factor 10~^ for the simple network and at least 10"^ for the 
complex network is need to recover the grain-free case. However, if metals are already 
depleted (by dust grains), we find that a depletion factor of 10~^ — 10~^ is just enough for 
both the simple and complex network to recover the grain-free result. This can also be seen 
in Fig. [5] at larger radius where all metals are effectively adsorbed. 

The differences between the networks are nontrivial. We cannot recover the results of 
the complex network by rescaling parameters in the simple one. Therefore a large array of 
species and reactions may be necessary to calculating the conductivity of protostellar disks. 

4.3. Two Grain Sizes 

There is no guarantee that the size distribution can be approximated by grains of a single 
size, because different kinds of reactions (e.g., adsorption and desorption) depend differently 
on grain size. Therefore, it is necessary to explore the chemistry in protostellar disks using 
more than one population of grains. On the other hand, the behaviors of grains with different 
sizes are almost independent. Neglecting grain coagulation, the direct interaction between 
the two grain populations is via charge exchange, which is slow (see §3.3. ip . As a result. 
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we expect to see that J2 ^i^^l (total grain surface area) or ^ NiQi roughly controls the free 
electron abundance, as we have found in the previous subsection. 

4.3.1. Chemistry 

Let Oi, /i be the size and mass fraction of the smaller grains, 02, /2 be the size and 
mass fraction of the bigger grains. Fixed total grain surface area means Niaf + A^2a2 
fi/di + 12/0.2 =const. Fixed A^'iai + A^202 means fi/a{ + /2/a2 =const. In Figure O we 
consider relatively small grains, ai = O.Ol/xm, 02 = 0.1/im. As in section 14.2. 1^ one sees 
that when grains are fully abundant, and when density is not too low, total grain surface 
area controls Xe. In the low density region, and when grains are substantially depleted, the 
controlling parameter of Xe lies in somewhere between ^ iVjOi and ^ Niof for the simple 
chemical network. For the complex network, Yl^i^i better controls Xe in these regimes. 
In Figure [3, we consider the two populations of grains to be larger, with ai = O.l/xm, 
02 = l.Oyum. We see the same trend as before. One notable difference is that for the bold 
solid lines in Fig. [Tb (i.e., complex model at fixed total surface area), there is a sharp increase 
as /2 10~^, or as /i —>■ 0, especially at low densities. This means that adding a population 
of smaller grains decreases Xe substantially. 

In all, the results above confirm the conclusions in §4.2. 1[ Small grains dominate the 
free electron abundance, especially in the complex network. 

4-3.2. Application to the Disk 

Figure [8] plots S(i(r) for various parameter choices. One can see that if the dust is fully 
abundant (/ = 10~^), the thickness of the active layer is appreciably less than lgcm~^. For 
the simple network, a reduction factor < 10~^ (/ < 10~^) is needed for the thickness of the 
active zone to reach 10 gcm~^. For the complex network, a reduction factor of 10~^ can raise 
Sactivc to 3gcm~^, whereas reduction by 10"^ (/ = 10~^) is still not sufficient to raise 
to lOgcm"^. Cosmic-rays hardly affect the ionization of the disk near 1 AU. The presence 
of the smallest dust grains (a = O.Ol/zm) demands a very high ionization rate of the order 

~ 10^^ to achieve = lOgcm"^ for standard values of the other parameters. This is 
excessively high for interstellar cosmic-rays but might perhaps be produced by nonthermal 
processes (reconnection?) within the corona of the disk or the protostar. 

We conclude that when we account for the smallest grains, the column density of the 
active zone is dramatically reduced. In order that the active column be large enough to 
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Fig. 6. — Like Fig. H] but for two grain populations, ai = 0.01/im and a2 = 0.1/im. Plots 
with bold lines have fixed total surface area /i/ai + /2/a2 =const, while plots with thin lines 
keep fi/af + /2/02 constant. Shown in the plots are Xg versus mass fraction of the larger 
grains, f2- (a). Bold lines: /i/ai + /2/a2 = O.lyum"^; thin lines: fi/al + f2/al = 1.0/im~^. 
(b). Bold lines: /i/ai + /2/«2 = 10~^/im~^; thin lines: /i/ai + /2/«2 = O.Olyum^^. Note that 
in the upper part of panel (b), thin and thick solid lines almost overlap. 
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Fig. 7. — Same as Fig. [6l but the two grain populations are replaced by ai = 0.1/xm, 
02 = 1.0/im. (a) Bold lines: /i/ai + /2/a2 = 0.01/im~^; thin lines: fi/ a{ + = 0.01/im~^. 
(b) Bold lines: /i/fli + /2/a2 = 10~^/im~^; thin lines: /i/fli + f2/0'2 = 10~^/im~^. 
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a) Simple chemical network 



10^ 



10 



" 10' 



10 



10" 



grain free ^ 




f=10"^ 




f=io-^ 

f=10-^ 


C„=10--s>>s 




— 




■ L^=5x10^°erg/s 

A 





10 

Radius (AU) 
b) Complex chemical network 



50 




Radius (AU) 



Fig. 8. — Like Figs. [2] & [5], but for two populations of grains with sizes oi = 0.01/im, 
a2 = 0.1/im. The mass ratio of the two populations is /1//2 = (02/0.1)^''^ = in rough 
accord with eq. (l23i) . We use / = /i + /2 to represent the total grain mass fraction. All 
parameters at standard values (Table [2]) except as indicated. 
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support vigorous accretion, the smallest grains must be severely suppressed. In the complex 
network, is more sensitive to the smallest grains, and more drastic depletion is needed. 
Our results show that O.Ol/xm grains should be almost completely removed, and 0.1/im ones 
should be reduced by a factor of IC^ compared to interstellar abundances. 



5. Molecular cooling of the active layer 

If accretion is driven by turbulent angular-momentum transport within the disk, then 



in steady state, the heat released per unit area from each side of the disk is (jPringlelll98ll ) 
Q ^ 3Mf2^/87r; the corresponding effective temperature is 



lbQM%^rll'^K. (30) 



The heat of wind-driven accretion would be less, because of the mechanical energy carried 
by the wind itself. At a normal interstellar gas-to-dust ratio, the heat of accretion would 
be radiated primarily by the dust. However, as we have seen, the dust abundance must 
be substantially lowered in the active layer in order that the conductivity of the layer be 
sufficient to couple it to the magnetic field. It is possible that the dust is so strongly depleted 
that the layer becomes optically thin. In that case, molecular lines might be seen in emission 
from the layer. 



In fact. lSalyk et al.l (120081 ) have observed H2O emission in the 10 — 20/im (with Spitzer- 
IRS) and 3 — 5/xm (with Keck-NIRSPEC) wavelength regions in two T Tauri systems with 
particularly high accretion rates of order 10~^^^Mq yr~^, DR Tau and AS 205. These authors 
infer that the emitting gas lies at ~ 1 AU from the central star in both systems and has a 
temperature ~ 10^ K. This is probably too hot to represent the active layer as a whole 
unless its optical depth is as small as r < lO^^ (j^^^g ^ r-^^^T^s if r < 1). Indeed, 
their analysis suggests a surface density < 0.1 gcm~^ for the emitting gas, assuming a cosmic 
abundance of oxygen in the form of H2O . So the emission may be coming from tenuous 
UV-heated gas well above the active layer. Nevertheless, these observations further motivate 
us to consider the possibility of significant molecular emission from active layers. 

In order that the molecular lines should dominate the cooling, however, the active layer 
must be more than modestly optically thin to dust, that is to say, it must be that Xd <C 1. 
The reason is that, in contrast to the situation in planetary or stellar atmospheres, the 
molecular lines are significantly narrower than their separation, so that they cover only a 
small fraction of the infrared spectrum. In other words, the emissivity of the active layer 
due to molecular lines is expected to be small. The rest of the present section is devoted to 
quantifying this statement via a representative calculation. 
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Since we are not attempting to fit actual data but only to indicate what might be 
expected to be emitted by a theorist's notional active layer, we consider H2O only. This 
is probably the most important molecular coolant because it is expected to be abundant, 
being composed of two of the most abundant elements and being strongly thermodynamically 
stable under these physical conditions, and because it has a much richer rotational spectrum 
than the common linear molecules CO and CO2. For similar reasons, ammonia (NH3) may 
be almost as important as water as a coolant. Because of the terrestrial significance of water 



vapor, the molecular spectrum of H2O is pa rticularly we 



on the extensive line and energy-level lists of (IBarber et al 



1 stud ied. We have relied mainly 



20061 . hereafter BT), but we have 



also consulted the JPL Molecular Spectroscopy database! and obtained consistent results 
for the molecular emissivity at Tgas = 300 K; however, the JPL database omits lines with 
wavelengths < lO/xm, whereas ~ 25% of a blackbody emits shortward of that cutoff even at 
300 K, and of course a larger fraction at higher temperatures. 

For simplicity, we represent the active layer by an isothermal slab in LTE at gas tem- 
perature T and define its emissivity by 



00 

^ / (l-e--)7r5.(T)rfz/, 



(31) 







where t^, is the frequency-dependent optical depth due to water alone, which in turn is related 
to the water column N (in molecules cm~^), the line "intensities" Ik (in cm molecule"^), 
and the line broadening function 0/c(Az/) (in cm) by 

T, = Nj2hMi'-''k). (32) 

k 

Following the molecular spectroscopists' convention, the frequency u is measured in wavenum- 
bers (cm~^) rather t han Hertz. We take i V/Sa = 3.6 x 10^° molecules g~^, which corresponds 



to the abundances of I Anders fc Gre vessel (119891 ) if all of the oxygen is in water. The relation- 



ship between the line intensity / and the Einstein coefficient Aij is temperature dependent 
and is given for LTE by eq. (3) of BT. 

In the limit of low columns where all lines are unsaturated (tj, ^ 1) eq. (13T!) would 
reduce to 

^- = ^E^'^^-(^)' (33) 

k 

This becomes 2.17 x 10~^°A^ at 300 K, corresponding to a Planck-mean opacity npi ^ 
7.8cm^g~^. Since the emissivity cannot exceed unity, it is clear that the important lines 



^spec . jpl . nasa . gov 
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must saturate at Sa > Kp/ ~ 0.13 gcm~^, and in fact the stronger lines saturate at even 
lower columns because under the conditions of interest, the lines are very narrow. Thus, for 
an active layer of surface density ~ lOgcm"^, the actual emissivity defined by eq. fl^ will 
be sensitive to the broadening prescription, which therefore merits some discussion. 
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Fig. 9. — Frequency- dependent absorption cross-section per water molecule over a small 
range of wavenumbers u = X^^ at T = 300 K and P = 1 yubar. Red curve: Thermal Doppler 
broadening only. Black: Doppler plus collisional broadening. Horizontal dashed line shows 
level above which lines would be saturated for a total mass column of lOgcm"^ assuming 



a solar abundance of oxygen in H2O, i.e. 3.6 x 10^*^ molecules g 
two overlapping lines; but note that one is relatively weak. 



Inset shows a closeup of 



The relative sizes of the natural width Az/q, the collisional width Auc, the thermal 
Doppler width Aud, and the turbulent width Aut (all of these are defined as half widths 
at half maximum) depend upon the temperature, pressure, and turbulent intensity. The 
temperatures of interest to us are within a factor of a few of that in the Earth's atmosphere, 
T ~ 300 X 10^°'^ K, but the pressures are much lower. If the base of the active layer 
lies at a height z^, <ti r above the midplane, then the pressure there is Pa ~ ^a^'^Za] if 
Zs, ~ 3[ksT / m{H2)Y^'^^~^ , i.e. three times the gaussian scale height of the disk (note that T 
here should be the temperature near the midplane, which may be somewhat lower than in 
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the active layer), then 

in other words, we are interested in pressures of order one microbar or less (10^ dyncm"^ = 
Ibar). The natural width 

^^^0 = ^^^^/' (35) 
/ 

is < 3 Hz (i.e., < 10~^°cm~^ in wavenumbers) for all upper levels i tabulated by BT whose 
energies are less than (3000K)A;b above the ground state. This is completely negligible 
compared to the other causes of line broadening; in particular, it is much smaller than the 
collisional width, which justifies our assumption of LTE. Neither BT nor the JPL database 
give collisional widths. These are difficult to calculate precisely and are well-measured for 
only a minority of transitions. For our purposes, it will be enough to have a rough estimate, 
so for all lines we use 



comparison. 



based on th e value quoted by iTownes fc Schawlowl (119551 ) for NH3 (not H2O ) in H2. For 



Giesen et all (11992h find Ai/, = 0.086 ± 0.002 cm"^ bar'^ for H2O in N2 at 
296 K. Thus, even at P ~ 1 /ibar, Ai/^ ^ IO^Az/q for the important transitions. So we have 
neglected the natural widths. The thermal Doppler width of a line with central frequency 
is 

1/2 . . 1/2 



Aud = V2 In 2 , , Uk ^ 1.46 x 10"^ cm"^ . (37) 

V^(H20)c2y Vsook; VlOOOcm-V ^ ^ 

This is typically two orders of magnitude smaller than the collisional width under terrestrial 
conditions, but at the much lower pressures of an active layer, Doppler broadening should 
dominate by a large factor: Auo ~ lO^Az/c. Nevertheless, the collisional broadening is 
not entirely negligible for us because it produces an approximately Lorentzian profile whose 
wings exceed those of the thermal Maxwellian far from the line center, and which can yield 
significant emission from the most strongly saturated lines. Finally, while the turbulent 
width is uncertain because of uncertainties in the strength and nature of the turbulence 
itself, for MRI-driven accretion we expect that the turbulent Mach number should be < 
0.5[q;A;b^/'^(H2)]^^^- This is less than the thermal width unless the viscosity parameter 
a > 0.5, which seems unlikely in these rather resistive circumstances. On these grounds, 
we have neglected the turbulent width. But it should be borne in mind that the turbulent 
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velocity fluctuations might be significantly nongaussian if the turbulence is intermittent, in 
which case the wings of the turbulent velocity profile might contribute importantly to the 
emission from the more saturated lines. Figure M demonstrates the effects of Doppler and 
collisional broadening on the frequency-dependent molecular absorption cross section, a,^; if it 
were due entirely to water, the corresponding opacity would be k^, = a^N^X (B.2O) / m(B.20) , 
where X(H20) is the abundance of water by mass, and m(H20) ~ ISrUp is the mass per 
molecule. At the very low pressures (by comparison with atmospheric or stellar conditions), 
relatively low temperatures, and modest columns relevant here, the strong and saturated 
lines rarely overlap. The Doppler-broadened spectrum {red online) shown in the figure was 
produced with the code spectra-bt2 of BT, and then convolved with a Lorentzian to 
simulate collisional broadening. 
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Fig. 10. — Cumulative emissivity of dust-free isothermal slabs due to water. Solid curves: 
E = lOgcm-2; dotted: Igcm'^. Top to bottom at right: T = 600,450,300, 150K. 

Figure [TO] shows the cumulative emissivity defined by the incomplete integral corre- 
sponding to eq. fl3Tl) . 

V 

e(<z.) = -^ 1 {\-e-^^')BAT)dv' , (38) 


for isothermal slabs representative of notional active layers. Evidently, the molecular emis- 
sivity of active layers is expected to be <^ 1, but perhaps not negligible compared to dust 
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if the latter is as strongly depleted as good magnetic coupling requires. The sublinear de- 
pendence on column density demonstrates that the emissivity is dominated by saturated 
lines. 

We ha ve found that eg. (|38|) can not be calculated reliably using the frequency-dependent 
opacities of Sharp fc Burrowa ( 2007 ) because these were intended for use in stellar and plane- 
tary atmospheres where collisional broadening is far stronger than in disks: the saturation of 
the stronger lines is therefore underestimated, so that the resulting emissivities are typically 
several times larger than those shown in F ig. [TUl We have not tried to use the even more 
recent tabulations of iFreedman et al.l (120081 ) but would expect similar difficulties, since those 
opacities were computed for minimum pressures of 300 /xbar. Rather than construct our own 
tables of k^, with the required resolution (finer than ~ 10~^), we have exploited the 

fact that the important lines are well separated to replace eq. (1351) with a hne-by-line sum 
that properly accounts for the broadening and saturation of the individual lines: 



e(< z/) ^ J] 5,,(r) j [l- exp [-iV40(z/' - z/,)]} dv' 



(39) 



lines k 



It is easily seen that eq. (13911 strictly overestimates eq. (l38il to the extent that lines overlap. 
But direct comparisons using limited wavenumber ranges indicate that the relative error of 
the line-by-hne sum ( l39l) is small (<^ 10^^) at our pressures and temperatures. 



5.1. Implication of conductivity constraints for molecular lines 

Figure [TT] summarizes our constraints on the dust in the plane of dust abundance (/) 
versus grain radius (a), for single-sized grains, or maximum grain radius (ctmax) for power-law 
size distributions, n{a)da oc a~^'^da. The conductivity constraint is based on the complex 
chemical network, with single species calculations. Parameters that lead > lOg cm~^ are 
considered as allowed. For the power-law size distribution, we have not tried to reproduce 
the complexities seen in Figures [31H1 The calculations show that the dependencies of the 
free-electron abundance, Xg, and of the column density of the active layer, Sa, are not 
accurately given by power laws in the dust-to-gas ratio (/) or grain radius (a). However, 
Figures [5], [6], & [7] do suggest that at gas densities relevant to the base of the active layer 
(p ~ 10~^^gcm~^), it is roughly the case that Xe varies with /'(a) and a in the combination 
J daf'{a)/aP with an exponent p somewhere between 1 and 2, where f'{a) oc a^N{a) is the 
differential mass fraction of grains with size a. Therefore, for the purposes of Fig. [TTl we 
have taken J^"™'"' daf'{a)/aP as the controlling parameter for the conductivity, with p = 3/2. 
To account for the power-law models with the MRN size distribution, we use a single-size 
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grain size[p.m] 

Fig. 11. — Constraints on grain size distribution and abundance in the MMSN at 1 AU. 
Abscissa is grain radius qq for single-size models, N{a) oc 6{a — ao), or Omax for MRN power- 
law size distribution ( 123|) with amin = O.Ol/xm. Light solid line: Conductivity constraint for 
single-size grains, derived using the complex network with = lOg cm~^; allowed models 
lie below line, shaded in light grey. Heavy solid line: Corresponding constraint for MRN size 
distribution, using the effective grain mass fraction of the smallest grains in equation (HOj) . 
Allowed region is shaded in dark grey. Light long-dashed lines: Loci of optical depth Tdust = 1 
for single-sized grains in active layer with Sgas = lOgcm"^ and Tgas = 150, 300, 600 K from 
top to bottom at left. Heavy long-dashed lines: Corresponding loci for MRN models. Dash- 
dotted lines: Corresponding loci for equal frequency-averaged dust and H2O emissivities, 
assuming MRN dust size distribution. Asterisk: Typical ISM values. 
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jrain model at the minimum grain size amin with an effective mass fraction 



( ^3/2 fl^max —2 J 

f( \ _ '-"'"''^^ Jamiu ^ Qmin + y/aminamax , , . 

I a ^l^da /Oinax 

"rnin 

The minimum radius is fixed at Omin = O.Ol/im, which we beheve to be conservative: 
if grains as small as this do exist in the active layer, then even smaller ones are probably 
present — perhaps all the way down to the PAH regime a ~ lOA — and such very small grains 
may dominate the conductivity. However, we have not calculated any ionization networks 
for grains smaller than O.Ol/xm. Values of Omax larger than lO^/im = 1 mm are not shown in 
Fig. [H] because our own rough estimates suggest that larger grains would precipitate out of 
the active layer even in the presence of residual turbulence in the underlying "dead" zone 



at th e level a^?, ^ 10 \ a value that we ded uce from (grain-free) simulations (ITurner et al. 



20071 : lOishi et al.ll2007l : iTurner fc Sanoll2008h . 



Also shown in Figure [TT] are the loci at which an active layer of surface density Sa = 
10 g cm~^ would be marginally optically thick to to dust, and the much lower loci at which the 
emissivities of the layer due to H2O and dust would be comparable. The dust opacities were 



derive d from the thermally-averaged optical efficiency factors computed by iDraine &: Lee 



(119841 ): we assumed a mixture of 58% silicate and 42% graphite grains by mass. The Figure 
shows that the absorption opacity depends mainly on dust mass fraction provided Omax ^ 
3/im, because emission and absorption then occur mainly in the dipole regime. The r^^st = 1 
curves are also rather insensitive to the temperature, at least over the range we consider 
(150 — 600 K). The equal-emissivity curves are more sensitive, however, because of the 
temperature dependence of the molecular emissivity seen in Fig. [TOl 

The main conclusion that we draw from Figure [TT] is that the active layer should be 
optically thin to dust. Therefore molecular emission lines should stand out above the dust 
continuum if the layer is heated by dissipation of MRI turbulence. The detectability of the 
lines will be reduced by the Doppler broadening associated with the orbital motion of the 
gas, which is much larger than thermal and collisional broadening, unless the disk is observed 
face on or spatially resolved (e.g. with an interferometer). The upper solid line shows that 
if grains smaller than a ~ l/xm are entirely absent, then there does exist a regime in which 
Tdust > 1 despite the ionization being sufficient for good coupling. The lighter solid and 
dashed lines show that if grains smaller than a ~ l/zm are entirely absent, then there does 
exist a regime in which Tdust ~ 1 despite the ionization being sufficient for good coupling. In 
the more plausible situation described by the power-law models, where small grains occur but 
in reduced abundance compared to ISM values, due both to growth of the largest grains and 
to reduction (or precipitation) of the overall dust mass fraction, the conductivity constraint 
requires rjust ^ 1. It seems even to require that cooling of the active layer is dominated by 
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molecules rather than dust, unless grains smaller than ~ O.l/xm are essentially completely 
excluded. 



6. Discussion 

The preceding sections have shown that our understanding of the interplay among grains, 
magnetic fields, and thermal emissions in protostellar disks suffers from many uncertainties. 
Some of these uncertainties might be reduced in the foreseeable future by further research. 
Others seem likely to plague astronomers for quite some time. 



6.1. Uncertainties in the Conductivity Calculation 

Our criterion of the active zone is based on the choice of the critical magnetic Reynolds 
number Re^"* = 100. As discussed in ^ our estimate may be conservative, and our calcu- 
lation may provide the upper limit of the size of the active zone. 

We have adopted the minimum-mass solar nebular disk model. This model is idealized 
in two aspects. It has an empirical surface density profile designed to match the mass distri- 
bution in the solar system. In the MMSN, the mass density scales linearly with the surface 
density coefficient Sq, but the disk scale height is independent of Sq, and therefore so is the 
fiux of X-rays and cosmic-rays impinging on the disk. In fact, we have experimented with 
different Sq, and found that Sa is essentially independent of Sq, although at higher surface 
density, the active layer of the disk resides at higher altitude. The MMSN is also assumed 
to be vertically isothermal. In reality, in the presence of small grains, the surface layers 
should be warmer than the midp lane even in a disk passively heated by stellar irradiation 
Jchiang fc GoldreichI IiQQtI . Il999h . and warmer still if turbulent accretion is concentrated 



in active layers. Higher temperature reduces the recombination rate and thereby increases 
the electron abundance, but only weakly (oc T~^/^, see §4.1.11) . Therefore, we expect that 
the size of the active zone is not affected much by the disk temperature profile. Molecular 
emission, however, could be strongly affected 

Major uncertainties attend the ionization model. The ionization rate is proportional 
to the X-ray luminosity and also depends on the X-ray temperature of the bremsstrahlung 
spectrum. A stellar wind may shield cosmic rays from the disk. In the outer disk. X-ray 
ionization is inefficient when small dust grains are present, and the ionization rate critically 
depends on cosmic rays. These effects are well studied in this paper. A related uncertainty 
comes from metal abundances. Metal atoms are effective electron donors. In the absence 
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of grains, the electron abundance is directly correlated with metal abundance. Since our 
results show that grains must be severely suppressed in order for a reasonable amount of the 
disk to be active, under such conditions, metals may be important. Our standard parameter 
assumes a relatively high metal abundance {xm = 1-25 x 10~^). Our results also suggest that 
the abundance is unlikely to be much lower than this, since an even more severe constraint 
on the dust abundance would then be required. Moreover, in the outer disk where the 
disk temperature is low, essentially all the metals are adsorbed onto dust grains. Unless 
cosmic-rays can provide the ionization, the outer disk is unlikely to be very active. 

Grain properties are another major source of uncertainty, as studied in detail in this 
paper. We devote the next subsection discussing the grain size distribution. Here, we 
emphasize that the simple and complex reaction networks show different response to grain 
abundance. Namely, in the complete absence of grains, or at the full ISM dust-to-mass ratio 
(/ = 0.01), the complex network produces more free electrons. However, for intermediate 
dust abundances, the simple network gives a larger Xe, sometimes eve n 10 times larger. Th e 



Oppenheimer & Dalgarno model, though popular for its simplicity (eg jTurner fc Sandl2008l ) 
should therefore be used with caution when grains are assumed. 



6.2. Grain size distribution 



The basic theme of this paper is to compare constraints on the grain abundance from 
the requirement of adequate magnetic coupling to those provided by infrared observations 
of thermal disk emission. For a given total mass / in grains per unit mass of gas, however, 
these two constraints depend differently on the sizes of grains, to the point that no useful 
comparison can be made if the size distribution is regarded as a completely free function. 
Thus, it is important to consider what constraints are available for this distribution. To 
bound the discussion, we take the distribution to be a truncated power law similar to the 
standard MRN form fl23l) but possibly with different lower and upper cutoffs (amin and Omax) 
and with a different exponent p, i.e. N{a) oc a^*'. 

CoUisional cascades involving macroscopic bodies such as asteroids or Kuijper-Belt ob- 
jects result in powerlaws with p ^ 3.5 when the cohesive strength is independent of size; 
the slope becomes steeper (i.e. larger p) if large bodies a re weaker than small ones, and 
shallower in the oppo s ite ca se (jO'Brien fc GreenbergI |2005| . and references therein). In line 
with this, iJones et al.l (119961 ) have proposed that the MSN distribution results from collisions 
between grains overrun by shocks. In their model, the grains are taken to have strengths 
similar to those terrestrial rocks and minerals, so that the minimum relative velocity between 
grains required for fragmentation is lkms~^. This is much larger than the likely collisional 
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velocities of submicron grains in protostellar disks. 



Many authors have c onsidered the colhsiona l agglo meration of grains in protostellar 
disks. In a recent study, iDuUemond &: DominikI (120051 ) demonstrate that if the sticking 
probability is taken to be high and fragmentation is ignored, then almost all grains smaller 
than a ~ 100 /im would disappear in less than 10^ yr, causing the disks to become optically 
thin. As they note, this contradicts the observation that many or most T Tauri disks remain 
optically thick, as judged by their spectral energy distributions at 10 — 100 /im, to ages of 
at least 10^ yr; the suggestion is that small grains must be replenished by fragmentation 
processes. 



In a recent review, iNatta et al.l (120071 ) discuss the observational evidence for grain 
growth and grain processing in protostellar disks. Millimeter observations indicate the pres- 
ence of millimeter or even centimeter-sized grains in many disks. On the other hand, infrared 
silicate features — especially crystalhne features near lO/zm — point to the persistence of mi- 
cron or submicron-sized grains even in the oldest disks. PAH features are also sometimes 
seen. But there is no clear evolutionary trend in the grain size distribution with disk (or 
rather stellar) age, which points again to active replenishment of the small-grain population. 



6.3. Turbulent Mixing 



We have calculated the free-electron abundance and conductivity as a function of local 
parameters: density, temperature, grain abundance, and ionization rate, etc. Since X-rays 
and cosmic rays penetrate only a limited column, the resulting electron abundance has a 
large vertical gradient. Any turbulent diffusion would mix electrons and ions vertically. The 
consequences for the width of the active layer depend upon the ratio of the timescales for 
mixing and recombination. If accretion is driven by turbulence and the turbulent transport 
coefficient for mass is comparable to that for momentum (i.e. to the turbulent viscosity), 
then absent dust grains, the diffusion may be faster than recombina tion; as a result, the 
dead zone may be reduced or even eliminated (jllgner fc Nelson! l2006bl ). 



Recently, [Turner et al.l (120071 ) performed 3D MHD simulations of protostellar disks with 
vertical structure, Ohmic resistivity, and time-dependent chemical evolution based on the 
simple (Oppenheimer-Dalgarno) reaction network without grains, allowing the reactants 
to be passively adv ected. Turbulen t mixi ng led to a weak coupling of the dead zone to 
the magnetic field. Illgner &: NelsonI (120081 ) found in similar simula tions that the effect o f 
turbulent mixing critically depend on gas-phase metal abundance. [Turner fc Sand (120081 ) 
included a uniform mass fraction of l/im dust grains. Despite rapid recombination on grain 



-44- 



surfaces, the radial magnetic field generated in the active layer diffused toward the disk 
mid-plane, causing some accretion there and rendering it "undead" . 

The simulations cited above adopted microphysical param eters very favorable to MR L 



For example, in the absence of grains, as in [Turner et al.l (120071 ) and lllgner fc Nelson! (120081 ). 



Turner & Sand ( 


2008) 


seen. Ilener & Nelson 



(l2008l ) chose a rather high protostellar X-ray luminosity (lO^^ergs"^) and temperature 
(5keV) . While all of these simulations would be described as "vertically stratified" in 
the parlance of MRI shearing-box work because they included vertical gravitational accel- 
erations, buoyancy forces were excluded by the use of isothermal equations of state. True 
stratification, as expected in real disks, may inhibit vertical mixing. 

Most importantly, the relevant component of the turbulent stress, T^'t'^ was found to 
be at best constant with height. Since, in steady state accretion, the accretion rate is 
proportional to the areal integral of this stress [see eq. ([E])], the "undead" zones contribute 
to M only in proportion to their geometrical thickness, which is comparable to that of the 
active layers, rather than in proportion to their mass column. 



7. Summary and Conclusion 



We have studied constraints on small dust grains in protostellar disks posed by the 
requirement that the electrical conductivity be sufficient to support magnetically driven 
accretion at observed rates. We have considered the implications of these constraints for the 
production of molecular emission lines, especially at radii ~ 0.1 — 10 AU where accretion is 
likely to to be confined near the surfaces of the disk in "active" layers comprising a fraction of 
the total column density. We adopt the minimum-mass solar nebular (MMSN) and assume 
protostellar X-rays and interstellar cosmic-rays t o be the main sources of i o nizati on. We 



compare a simple chemical reaction modeled on lOppenheimer fc Da. 



garnol (1974 ) with 



complex reaction network extracted from the latest UMIST database (jWoodall et al.ll2007l ). 
Reactions involving du st grains are modeled according to the general scheme laid out by 
Ilgner fc NelsonI (j2006al . hereafter IN06a), with some improvements, including the use of up 
to two grain sizes in order to mimic an MRN grain size distribution. The conductivity is 
heavily dominated by the smaller grains. We discuss in detail the dependence of free electron 
abundance (xe) on density, temperature, ionization rate, gas-phase metal abundance, and 
especially on grain size and grain abundance. 
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Separately, we estimate that electron fractions Xe ^ 10~^^-10~^° are required at 1 AU 
to explain accretion at 10~'^ MqJt~^ , which is typical of the more active systems; minimum 
magnetic field strengths range from a few tenths of a Gauss if accretion is driven by winds, 
to several gauss if it is powered by magnetorotational turbulence. It is difficult to be very 
precise about these requirements because of the complexities of winds and turbulence in the 
presence of vertical stratification and tensorial conductivities, but the numbers just given 
are probably conservative. We define the active layer as the column density at which 
R^M > 100, where Rcm = Cs/flr] and r] is the Ohmic diffusivity. Smaller values of Rcm are 
probably not realistic because, at the field strengths we estimate. Hall and ambipolar drift 
contribute to slippage between the neutrals and the field. By this definition, it is difficult to 
achieve > lOgcm"^ at 1 AU even if the small grains are strongly suppressed. 

In more detail, our conclusions regarding the conductivity are as follows: 



The free electron abundance Xe depends more strongly on density (p) and ionization 
rate {() than on temperature. Even a small amo unt of dust suppresses a tomic ions 



such as Mg"*" and Fe^ in the gas phase, as found by lllgner fc Nelson! (j2006af ): however, 
we find that this is due to recombination on grain surfaces rather than adsoprtion, 
except at temperatures (< lOO-ft' for Mg). 

In the absence of grains, the complex chemical reaction network predicts a slightly 
higher Xg than the simple Oppenheimer-Dalgarno model, rather than the reverse as 
IN06a found. This appears to be caused by updates to the reaction rates in the 
UMIST database. Th ese changes are within the typical uncertainties of those rates 



fIVasvunin et al.ll2008f ) 



The electron abundance depends sensitively on grains, but this dependence cannot be 
accurately characterized by a simple combination of grain size and grain abundance, 
especially in the complex network. Roughly however, the controlling parameter lies 
somewhere between the total grain surface area, ^a^A^(a), and the grain abundance 
weighted by linear size, ^aiV(a). 

Standard cosmic-ray ionization rates (cr = 10~^^-10~^^s~^ render the entire disk col- 
umn active beyond 10 AU even in the presence of sub-micron grains. X-rays contribute 
to Xe only near the disk surface, within columns S < lOgcm"^. The effect of X-ray 
scattering on the ionization rate is small at 1 AU, but becomes important farther out. 

In the inner disk around 1 AU, rather extreme parameters are required to achieve active 
layers as large as = lOgcm"^ if submicron grains are present. For example, if all 
grains have size a = 0.1/xm, then the dust-to-gas ratio must be reduced to / < 10"^ 
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for standard X-ray and cosmic-ray parameters. If (cr were increased to lO^^^s"^ — 
well above ISM values but perhaps a proxy for nonthermal processes within the disk 
or stellar corona — then / < 10~^ could be tolerated, still two orders of magnitude 
below the ISM abun dance but compatible with some models of the infrared spectral 



energy distribution (ID'Alessio et al J 120061). If the grains extend to a < O.Ol/xm with 



a standard N{a) oc a size distribution, then / < 10 ^ even if the maximum grain 
size grows to ~ 1 mm with our standard ionization sources. 

The last two points above lead us to suspect that either an important source of ionization 
has been overlooked, or else something other than MRI turbulence drives accretion, at least 
near 1 AU. However, if we suppress these suspicions, then we conclude that the active layer 
at 1 AU should be optically thin to dust in the Planck average. Just how thin depends 
upon the minimum grain size, since the dust opacity depends mainly on total dust mass 
fraction / as long as ctmax is smaller than a few microns. We have estimated the emissivity 
of the active layer due to water lines, taking into account that the important lines are 
strongly saturated for > 1 g cm^^ because collisional broadening is almost negligible at 
relevant pressures (< lyubar). For plausible grain properties, the heat dissipated by MRI 
turbulence may be radiated primarily in molecular lines, and the gas temperature may be 
significantly higher than the effective temperature because the total emissivity of the active 
layer due to dust and molecules is small. This offers the exciting prospect that molecular 
lines from a c tive l ayers may be observable, and indeed may already have been observed by 



Salyk et al.l (120081 ). However, quantitative diagnoses of active layers using such lines will 
require theoretical modeling of the coupled dynamical and thermal evolution of the MRI 
turbulence. 



We thank Bruce Draine for patiently educating us about grains and their interactions, 
M. Ilgner for advice concerning the integration of large reaction networks, and S. Cazaux for 
discussions of hopping rates on grain surfaces. We also thank Mark Wardle, and our referee, 
Neal Turner, for stressing the importance of X-ray scattering on disk ionization. This work 
was supported in part by NSF award PHY-0821899 "Center for Magnetic Self-Organization 
in Laboratory and Astrophysical Plasmas." 



A. Comments on X-ray Ionization 



The X-ray ionization rate is a crucial factor on disk conductivity. In our original ver- 
si on of this paper , we c alculated the X-ray using the formula given by equations (2)- (4) 
of iFromang et al.l ((2002), where the X-ray photons was assumed to be attenuated only by 
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absorption. However, at the energies of interest, the Compton cross section for photons is 
comparable to photoionization cross section. Since X-ray photons incident the protostellar 
disks obhquely, scattered photons can penetrate deeper inward, resulting in higher ioniza- 
tion rate toward disk middle plane. As supplemental information to section 3.1, we provide 
a comparison of ionization rate between pure absorptio n model and the model including 
the effect of Compton scattering ( Igea fc Glassgoldlll999l ). as shown in Fig. [121 This figure 
clarifies that neglecting Compton scattering would significantly underestimate the ionization 
rate at large radii, especially towards the disk midplane, but is OK at small radii around 

1 AU. 



-48- 





Fig. 12. — X-ray ionization rate as a function of disk column density, at 1 AU and 10 AU 



respectively. The central X-ray source is assumed to have X-ray luminosity Lx = 10 erg 
s^\ and X-ray temperatu re Ty = 3keV (u p per p anel) and Tx = 5keV (bottom panel). 
Circles: values taken from llgea fc Glassgoldl (119991 ). Asterisks: ionization rate with X-ray 
absorption only. Solid line: our fitted ionization rate for lAU. Dash ed line: rescaling of th e 
fitted curve to lOAU. Pink curves are the fitting formula adopted by lTurner fc Sand (120081 ). 
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ParameterMeaning Standard Value Range 

Stellar Mass IM© Fixed 

Eq Disk Surface Density at lAU 1700g cm~^ Fixed 

To Disk Temperature at lAU 280K Fixed 

Rbm Critical Magnetic Reynolds Number 100 >100 

Lx ProtoStellar X-ray Luminosity 0.5 x lO^^erg lO^^^^^erg s"^ 

Tx ProtoStar X-ray Temperature 3keV 1 — 5keV 

Co Cosmic-ray Ionization Rate Os~^ < 10~^^s~^ 

01 Radius of Small Grains 0.01/xm > 0.005/im 

02 Radius of Big Grains O.l/xm < 3/im 

fi Mass Fraction of Small Grain Particles 0.0076 — 0.01 

/a Mass Fraction of Big Grain Particles 0.0024 - 0.01 

Mass Density of Grains 3g cm~^ Fixed 

77 H2 Formation Efficiency 10~^ < 1 

xug Abundance of Mg 1.0 x 10"^ < 10"'^ 

xpe Abundance of Fe 2.5 x lO"'^ < 10"^ 



Table 2: Model parameters. 
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